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1.  INTRODUCTION 


This  report  documents  a  computational  aerothermal  analysis  of  the  XM797  gas  generator 
nose  cap.  The  XM797  nose  cap  is  a  critical  component  in  the  design  of  a  120-mm  finned 
kinetic  energy  (KE)  training  round.  The  concept  is  illustrated  in  Figure  1,  taken  from  Kotar 
and  Quevedo  (1983).  At  launch,  the  nose  cap  serves  as  a  forward  clamp,  holding  together 
multiple  sections  of  the  long-rod  projectile.  A  propellant-filled  chamber  within  the  nose  cap 
is  designed  so  that  aerodynamic  heating  ignites  the  propellant  at  some  distance  down  range, 
separating  all  the  major  projectile  components.  The  training  round  then  breaks  apart  into 
high  drag  pieces  that  quickly  decelerate  and  impact  the  ground.  Ballistic  match  is  desired 
between  the  fielded  and  training  rounds  to  a  range  of  3.5  km,  with  a  maximum  range  less 
than  12  km. 

The  objective  is  to  provide  the  propellant  temperature  response,  which  is  necessary  to 
design  the  chamber  geometry  for  training  applications.  A  model  for  the  aerodynamic  heat 
transfer  is  constructed  using  a  numerical  boundary-layer  technique.  This  model  provides  the 
necessary  boundary  conditions  for  a  numerical  heat  conduction  analysis,  which  computes  the 
temperature  response  throughout  the  nose  cap.  The  temperature  response  at  the  location  of 
propellant  ignition  is  the  key  parameter  of  interest.  Once  the  temperature  response  is  known, 
the  nose  cap  can  be  designed  for  the  propellant  (and  corresponding  ignition  temperature)  of 
choice. 

A  105-mm  version  of  the  XM797  nose  cap  was  previously  studied  (Abbett,  Duiven,  Laub, 
and  Beck  1981;  Kotar  and  Quevedo  1983;  Sturek,  Kayser,  and  Mylin  1983;  Sturek,  Kayser, 
Mylin,  and  Hudgins  1983).  For  validation,  analysis  for  the  105-mm  training  round  is  first 
presented  and  compared  with  the  results  from  Kotar  and  Quevedo  (1983).  Next,  results 
are  presented  for  the  current  120-mm  training  round  projectile.  The  modeling  approaches 
for  the  heat  transfer  and  heat  conduction  analyses  are  described,  and  the  relevant  modeling 
assumptions  axe  discussed. 


2.  CONFIGURATIONS  AND  FLIGHT  CONDITIONS 

Figure  2  shows  the  105-mm  version  of  the  XM797  nose  cap.  This  figure  was  reproduced 
from  the  report  by  Kotar  and  Quevedo  (1983),  in  which  it  is  referred  to  as  the  “revised 
design.”  The  nosetip  material  is  given  as  18-8ph  stainless  steel,  whose  thermal  properties  are 
assumed  to  be  the  same  as  304  stainless  steel,  shown  in  Table  1.  The  propellant  chamber 
is  a  right-circular  cylinder  situated  along  the  projectile  axis  and  is  sealed  with  a  Teflon® 
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cup.  The  nose  cap  is  secured  to  the  multiple-piece  rod  by  a  steel  piston/aluminum  shear  pin 
mechanism  designed  to  fail  upon  propellant  ignition.  The  in-flight  velocity  profile  (leading 
to  propellant  ignition)  is  assumed  to  be  that  of  the  105-mm  APFSDS-T  (armor-piercing,  fin- 
stabilized,  discarding  sabot  with  tracer)  M833  projectile,  as  given  in  Department  of  the  Army 
Firing  Table  FT  105-A-3.  The  launch  velocity  is  1.49  km/s,  and  the  velocity  retardation 
is  approximately  54.0  m/s/km.  For  the  heat  conduction  computations,  the  free  stream  and 
initial  wall  temperatures  are  assumed  to  be  324  K. 

Table  1.  Thermal  Properties  for  AIAI-SAE  304  Grade  Stainless  Steel 


Density,  p  (kg/m^) 

7850.0 

Conductivity,  k  (W/m/K) 

10. 

Specific  Heat  Cp  (J/kg/K) 

43.7 

Diffusivity,  (i^^/s) 

2.9x10-® 

Figure  3  shows  the  120-mm  version  of  the  XM797  nose  cap.  The  basic  components 
of  the  120-mm  nose  cap  are  the  same  as  those  of  the  105-mm  nose  cap.  The  physical 
dimensions  correspond  to  those  given  by  the  drawing  entitled  “Nose  cap  TM-2376,”  provided 
by  Mr.  Brian  C.  Wong,  U.S.  Army  Armament  Research,  Development,  &  Engineering 
Center,  Picatinny  Arsenal,  New  Jersey.  The  materials  are  assumed  to  be  the  same  as  those 
in  the  105-mm  version.  The  in-flight  velocity  profile  is  taken  to  be  that  of  the  120-mm 
APFSDS-T  M827  and  APFSDS-T  DM13  projectiles,  whose  profiles  are  almost  identical,  as 
given  in  Department  of  the  Army  Firing  Table  FT  120-D-l.  The  launch  velocity  is  1.65 
km/s,  and  the  velocity  retardation  is  approximately  74.0  m/s/km.  For  the  heat  conduction 
computations,  the  free  stream  and  initial  wall  temperatures  are  assumed  to  be  equal.  Three 
different  cases  are  considered:  272  K,  294  K,  and  324  K,  referred  to  as  the  “cold,”  “standard,” 
and  “hot”  cases,  respectively. 

The  angle  of  attack  is  prescribed  as  0°  for  all  cases.  Atmospheric  sea-level  free  stream 
pressure  and  density  are  assumed.  The  state  of  the  boundary  layer  (i.e.,  laminar,  turbulent, 
and/or  transitional)  is  discussed  individually  for  the  105-mm  and  120-mm  versions  in  the 
results  section. 


3.  MODELING  APPROACH 

3.1  Surface  Heat  Transfer  Model.  Surface  heat  transfer  characteristics  were  com¬ 
puted  using  an  axisymmetric  boundary-layer  numerical  technique  known  as  the  ASCC 
(ABRES  Shape  Change  Code)  reported  by  Suchsland  (1980).  The  ASCC  code  uses  a  mod- 
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ified  Newtonian  surface  pressure  distribution  and  inviscid  flow  field  approximations  to  gen¬ 
erate  boundary-layer  edge  flow  conditions.  The  ASCC  code  was  used  in  this  study  because 
it  generates  adiabatic  wall  temperatures  and  heat  transfer  coefficients  in  1  or  2  minutes  of 
central  processing  unit  (CPU)  time  on  a  typical  workstation  or  mini-computer,  as  opposed 
to  Navier-Stokes  methods,  which  usually  take  more  CPU  time. 


Previous  work  at  the  U.S.  Army  Research  Laboratory^  (ARL)  included  comparison  of 
ASCC  and  Navier-Stokes  results  for  KE  projectile  nosetips  and  published  experimental  mea¬ 
surements  (Guidos  and  Sturek  1993).  The  computed  laminar  heat  transfer  rates  compared 
well,  but  the  turbulent  results  showed  that  the  ASCC  code  underpredicted  the  heat  transfer 
rates  by  about  40%  at  1.7  km/s,  compared  to  the  Navier-Stokes  results.  For  this  reason,  the 
turbulent  heat  transfer  rates  computed  by  the  ASCC  code  for  this  study  were  increased  by 
40%  before  being  used  in  the  heat  conduction  analysis. 


The  formulation  of  the  heat  transfer  model  follows  the  procedure  outlined  by  Guidos  and 
Weinacht  (1993).  Newton’s  Law  of  Cooling  is  applied  to  form  a  distribution  of  heat  transfer 
coeflRcient,  h,  at  various  velocities  of  interest,  i.e.. 


h  = 


Ta^  -  n. 


(1) 


in  which  Taiu  is  the  computed  adiabatic  wall  temperature,  q  is  the  computed  heat  transfer 
rate,  and  is  wall  temperature  at  which  the  heat  transfer  rate  is  computed.  In  the  study 
done  by  Guidos  and  Weinacht  (1993),  Taw  was  held  constant  over  the  body  surface  at  each 
velocity,  since  the  nosetip  was  ignored  in  the  analysis.  In  the  present  study.  Taw  is  computed 
as  a  function  of  axial  distance  from  the  nosetip  (i.e.,  stagnation  point)  at  two  different 
velocities  bound  that  of  the  projectile  in  flight. 


The  adiabatic  wall  temperature,  Taw,  is 
defined  as 


represented  by  a  free  stream  recovery  factor,  r/. 


Taw  -  Too 
To  —  Too 


(2) 


in  which  Too  and  To  are  the  free  stream  and  total  temperatures,  respectively.  The  use  of  r/ 
serves  to  normalize  the  behavior  of  To^,  for  different  flight  velocities. 


For  the  105-mm  case,  the  heat  transfer  model  is  formed  using  results  at  velocities  1.49 
km/s  and  1.27  km/s;  for  the  120-mm  case,  the  heat  transfer  model  is  formed  using  results  at 
velocities  1.65  km/s  and  1.43  km/s.  In  both  cases,  the  ASCC  results  are  generated  assuming 
a  free  stream  temperature  of  294  K.  For  convenient  reference,  the  total  temperatures  for 
each  of  the  velocities  of  interest  is  calculated  from  compressible  flow  theory  (Zucrow  and 
Hoffman  1976)  and  is  given  in  Table  2. 

'  Formerly  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL) 
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Table  2.  Total  Temperature  for  Flight  Velocities  of  Interest,  Assuming  Too=294  K 


Velocity  (km/s) 

Mach  Number 

r„(K) 

1.65 

4.8 

1.49 

4.3 

1.43 

4.16 

1312. 

1.27 

3.7 

1099. 

The  heat  transfer  coefficients  are  held  constant  with  respect  to  wall  temperature  and 
free-stream  temperature.  In  the  study  of  Guidos  and  Weinacht  (1993),  it  was  found  that  h 
decreases  slightly  as  increases  toward  Taw  Thus,  the  assumption  of  constant  h  based  on 
atmospheric  wall  temperatures  will  contribute  to  a  slight  underprediction  of  the  projectile 
temperature. 


3.2  In-Depth  Heat  Conduction.  The  transient  in-depth  heat  conduction  analysis 
was  made  using  the  numerical  code  described  by  Dwyer  (1990);  Sturek,  Dwyer,  and  Ferry 
(1990);  Yam  (1991);  and  Sturek  (1993).  The  scheme  is  a  3-D,  time-accurate,  implicit,  finite 
volume  technique  cast  in  generalized  coordinates.  For  this  study,  an  axisymmetric  version  of 
the  code  was  used.  The  computations  were  performed  using  a  single  processor  on  a  Silicon 
Graphics  Inc.  Challenge  XL/R4400  computer,  with  each  case  requiring  approximately  1.5 
hours  of  CPU  time.  One  of  the  strengths  of  this  code  is  its  fast  execution,  because  of  its 
implicit  formulation.  Some  weaknesses  of  the  code  are  that  its  application  is  limited  to 
single-zone,  single-material,  constant  thermal  property  cases.  The  technique  has  recently 
been  validated  for  the  case  of  axisymmetric  convective  heating  of  an  infinite-length  rod,  for 
which  an  analytical  solution  (Schneider  1955)  is  available  (see  the  appendix). 

The  problem  setup  is  shown  in  Figure  4,  and  the  grid  is  shown  in  Figure  5.  The  cup  seal 
is  used  as  the  aft  grid  boundary  and  is  assumed  to  form  an  adiabatic  surface  perpendicular  to 
the  projectile  axis.  Table  3  shows  the  thermal  properties  of  Du  Pont  Teflon®,  obtained  from 
a  Du  Pont  Product  Information  Chart  (1995)  and  Modern  Plastics  Encyclopedia  (1991). 
The  diffusivity  is  two  orders  of  magnitude  less  than  that  of  stainless  steel,  demonstrating  the 
insulative  nature  of  the  seal.  The  adiabatic  boundary  formed  by  the  seal  is  then  assumed  to 
extend  radially  outward  to  the  projectile  surface. 

Conforming  to  the  single-zone,  single-material  restriction,  the  propellant  chamber  is  mod¬ 
eled  as  solid  stainless  steel.  Since  the  chamber  is  located  in  the  general  area  of  lowest  expected 
thermal  gradients,  this  simplification  is  acceptable  as  a  baseline  case.  But  the  thermal  behav¬ 
ior  of  the  propellant  is  worth  considering  in  more  detail.  Table  4  shows  approximate  thermal 
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Table  3.  Thermal  Properties  for  Du  Pont  Teflon® 


Density,  p  (kg/m^) 

1600. 

Conductivity,  k  (W/m/K) 

0.239 

Specific  Heat  Cp  (J/kg/K) 

1088. 

Diffusivity,  o=~  (n^^/s) 

1.4x10“'^ 

properties  of  several  propellants  following  Cohen,  Beyer,  Bilyk  and  Newberry  (1992)  and 
Ward  (1974).  The  properties  are  quite  similar  to  those  of  the  cup  seal,  demonstrating  that 
the  propellant  chamber  serves  as  a  good  insulator  for  the  rest  of  the  nose  cap  section.  It  may 
be  concluded  that  modeling  the  propellant  chamber  as  solid  stainless  steel  will  contribute  to 
a  slight  underprediction  of  the  nose  cap  temperature. 

Table  4,  Thermal  Properties  for  M9,  JA2,  M30,  and  X14  Propellants 


Density,  p  (kg/m^) 

1600. 

Conductivity,  k  (W/m/K) 

0.234 

Specific  Heat  Cp  (J/kg/K) 

1465. 

Diffusivity,  (j^^/s) 

1.0x10-'^ 

In  retrospect,  an  improved  prediction  could  probably  be  obtained  from  the  single-zone, 
single-material  approach  by  modeling  the  chamber  wall  as  adiabatic  with  respect  to  the 
heat  flow  in  the  stainless  steel.  The  penalty  would  be  the  need  to  construct  an  irregular 
computational  grid  coincident  with  the  chamber  wall,  rather  than  the  current  grid  which  is 
coincident  with  the  projectile  axis.  An  effort  is  under  way  to  use  a  newly  available  multiple 
zone,  multiple  material,  variable  thermal  properties  heat  conduction  code. 

4.  RESULTS 

4.1  105-mm  XM797  Nose  cap.  As  a  means  of  validating  the  current  modeling 

approach  for  a  similar  problem,  analysis  was  first  conducted  for  the  105-mm  version  of  the 
XM797  nose  cap,  for  which  previous  results  have  been  published  by  Kotar  and  Quevedo 
(1983).  Few  details  of  the  heat  transfer  and  heat  conduction  modeling,  which  was  conducted 
by  AVCO,  are  included  in  the  aforementioned  report.  However,  that  report  does  include 
some  additional  nosetip  transient  thermal  response  predictions  that  compare  favorably  with 
unreferenced  arc  jet  tests,  giving  a  measure  of  confidence  to  the  free  flight  results  presented 
there. 
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In  the  current  approach,  the  aerodynamic  heat  transfer  parameters  are  computed  using 
the  ASCC  code  before  computing  the  in-depth  heat  conduction.  Figure  6  shows  the  com¬ 
puted  recovery  factors  for  the  105-mm  XM797  nose  cap  at  1.49  km/s  (launch  velocity)  and 
1.27  km/s  (uninterrupted  velocity  at  4-km  range).  The  distribution  is  shown  as  a  function 
of  axial  location,  x,  normalized  by  a  reference  length,  /,  equal  to  25.4  mm  (1  inch).  The 
recovery  factor  as  computed  by  the  ASCC  code  is  the  same  for  laminar  and  turbulent  flows. 
Much  experimental  data  exists  showing  the  recovery  factor  to  differ  between  laminar  and 
turbulent  flows  by  about  5%  (see  Guidos  &  Weinacht  [1993]  for  a  more  detailed  discussion). 
In  the  present  study,  the  accuracy  of  analysis  is  not  compromised  because  the  same  recovery 
factor  distribution  is  used  in  both  the  formulation  and  application  of  the  heat  transfer  model. 

Figures  7  and  8  show  the  computed  105-mm  XM797  heat  transfer  coefficients  used  in 
the  present  heat  conduction  analysis.  The  fully  turbulent  and  fully  laminar  cases  are  shown. 
As  previously  mentioned,  the  turbulent  heat  transfer  coefficients  computed  by  the  ASCC 
code  were  scaled  by  a  factor  of  1.4  before  being  used  in  the  heat  conduction  analysis.  It  is 
assumed  that  the  heat  transfer  coefficient  and  adiabatic  wall  temperature  both  vary  linearly 
with  respect  to  velocity  during  the  trajectory. 

Figure  9  shows  the  computed  temperature  response  at  the  location  of  propellant  ignition 
within  the  nose  cap.  Propellant  ignition  is  assumed  to  occur  at  the  forward  outer  corner 
of  the  propellant  chamber.  The  forward  corner  of  the  propellant  chamber  is  located  19.94 
mm  (0.785  inches)  down  stream  from  the  stagnation  point  and  2.51  mm  (0.099  inches)  off 
axis,  as  shown  in  Figure  2.  The  present  results  are  shown  for  both  the  fully  laminar  and 
fully  turbulent  heating  cases.  The  computed  temperature  difference  between  the  laminar 
and  turbulent  cases  is  near  150  K  at  2  seconds  into  the  flight,  indicating  large  sensitivity  of 
propellant  ignition  to  the  state  of  the  boundary  layer.  Considering  that  a  typical  propellant 
ignition  temperature  is  500  K  or  more,  it  may  be  concluded  that  a  laminar  flow  situation 
will  not  heat  the  nose  cap  sufficiently  to  ignite  the  propellant,  making  a  boundary  layer  trip 
(i.e.  surface  knurling)  a  necessary  design  component. 

Figure  9  also  includes  the  results  from  Kotar  and  Quevedo  (1983),  which  differ  from 
the  present  turbulent  results  by,  at  most,  only  20  K  during  the  first  4  seconds  of  flight. 
The  comparison  suggests  that  the  previously  published  results  were  generated  for  a  fully 
turbulent  flow  case. 


4.2  120-mm  XM797  Nose  cap.  With  the  modeling  approach  validated  for  the  105- 

mm  version  of  the  XM797  nose  cap,  results  were  then  generated  for  the  120-mm  version. 
Figure  10  shows  the  computed  recovery  factors  for  the  120-mm  version  at  1.65  km/s  (launch 
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velocity)  and  1.43  km/s  (uninterrupted  velocity  at  3  km  range).  As  was  the  case  for  the 
105-mm  version,  the  recovery  factors  computed  by  the  ASCC  code  were  the  same  for  laminar 
and  turbulent  flows. 

Figures  11  and  12  show  the  computed  120-mm  XM797  heat  transfer  coefficients  used  in 
the  heat  conduction  analysis.  The  analysis  was  done  by  assuming  boundary  layer  transition 
to  occur  over  a  section  of  surface  knurling  (not  shown  in  Figure  3)  near  the  nosetip.  The 
heat  transfer  coefficient  in  the  transition  region,  htr,  was  generated  by  combining  the  fully 
laminar  and  fully  turbulent  distributions  using  the  expression 

htr  =  hi  +  Kt{ht  —  hi)  for  Xi  <  x  <  X2  (3) 

in  which  hi  and  ht  are  the  laminar  and  turbulent  heat  transfer  coefficients,  respectively,  at 
a  specified  x  location;  X\  (equal  to  2  mm)  and  X2  (equal  to  10  mm)  are  the  beginning  and 
end,  respectively,  of  transition;  and  is  the  transition  intermittency  factor.  This  factor  is 
identically  zero  for  laminar  flow  and  identically  one  for  fully  turbulent  flow.  The  aissumed 
form  of  Kt  in  the  transition  region  is  that  of  a  5th  order  polynomial,  as  previously  used  by 
Guidos,  Weinacht,  and  Dolling  (1990).  As  in  the  105-mm  case,  the  turbulent  heat  transfer 
coefficients  computed  by  the  ASCC  code  were  scaled  by  a  factor  of  1.4  before  being  used  in 
the  heat  conduction  analysis. 

Figure  13  shows  the  computed  temperature  response  at  the  location  of  propellant  ignition 
within  the  nose  cap.  Once  again,  propellant  ignition  is  assumed  to  occur  at  the  forward  outer 
corner  of  the  propellant  chamber.  The  forward  corner  of  the  propellant  chamber  is  located 
23.62  mm  (0.93  inches)  down  stream  from  the  stagnation  point  and  2.3  mm  (0.0905  inches) 
off  axis,  as  shown  in  Figure  3.  For  a  propellant  ignition  temperature  of  500  K,  ignition  is 
computed  to  occur  at  about  2.1  seconds  (3  km  range)  into  flight.  The  cold,  standard,  and 
hot  cases  are  shown  to  differ  by  a  total  of  30  K  or  less  for  flight  times  greater  than  2  seconds. 
This  translates  into  less  than  1/3  of  a  second  variation  in  propellant  ignition  onset. 

Figure  14  shows  computed  temperature  contours  within  the  nose  cap  at  several  times 
during  flight  for  the  standard  case.  The  orientation  of  the  contours  verifies  the  contention 
that  the  highest  chamber  temperature  (thus  ignition)  occurs  at  the  forward  corner  of  the 
chamber. 

Figure  15  presents  the  complete  temperature  distribution  at  various  times  in  a  form 
useful  for  possible  design  modifications.  The  temperature  distribution  is  shown  as  a  function 
of  propellant  chamber  depth,  djl.  The  chamber  depth  is  measured  from  the  base  of  the 
nose  cap  (refer  to  Figure  3)  and  is  currently  24.89  mm  (0.98  in).  Figure  14  depicts  the 
fact  that  the  forward  corner  of  a  deeper  propellant  chamber  is  closer  to  the  nosetip,  and 
accordingly  experiences  a  higher  heating  rate.  Noting  that  the  time  of  flight  at  3.5  km  range 
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is  approximately  2.3  seconds,  the  depth  of  the  chamber  can  be  varied  to  correspond  to  the 
selected  propellant’s  ignition  temperature.  Figure  15  shows  that  for  a  propellant  ignition 
temperature  of  500  K,  a  2.54  mm  (0.1  inch)  increase  in  chamber  depth  hastens  the  ignition 
time  by  slightly  less  than  1/2  second. 


5.  CONCLUSION 

A  computational  aerothermal  analysis  of  the  XM797  gas  generator  nose  cap  was  pre¬ 
sented.  The  temperature  response  of  the  nose  cap  is  a  critical  component  in  the  design  of 
105-mm  and  120-mm  finned  KE  training  round  projectiles,  which  use  the  XM797  concept. 
The  temperature  response  within  the  propellant  chamber  inside  the  nose  cap  is  the  key 
parameter  of  interest. 

Analysis  was  first  presented  for  a  105-mm  version  of  the  XM797  nose  cap  using  flight 
conditions  from  an  M833  projectile.  The  computed  temperature  response  at  the  forward 
corner  of  the  propellant  chamber  compared  favorably  with  previously  published  results, 
validating  the  present  modeling  approach. 

Analysis  for  the  120-mm  version  of  the  XM797  nose  cap  was  performed  using  flight 
conditions  from  a  DM13  projectile.  The  analysis  was  done  for  the  boundary-layer  transition 
case,  in  which  the  boundary  layer  is  tripped  by  a  section  of  surface  knurling  near  the  nosetip. 
For  a  propellant  ignition  temperature  of  500  K,  ignition  is  computed  to  occur  at  about  2.1 
seconds  (3  km  range)  into  the  flight.  Results  for  cold,  standard,  and  hot  wall  cases  show 
relatively  minor  differences  which  translate  into  less  than  1/3  second  variation  in  propellant 
ignition  onset.  For  the  current  design,  a  2.54-mm  (0.1-inch)  increase  in  chamber  depth 
hastens  the  ignition  time  by  approximately  1/2  second. 
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Figure  1.  XM797  Gets  Generator  Nosecap  Concept 


1.91 


Dimensions  in  inches 


Figure  2.  Schematic  of  105-mm  XM797  Nosecap 
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Projectile  Axis 


Propellant  Chamber 


Figure  4.  Heat  Conduction  Problem  Setup 


Figure  5.  Single- Zone,  Single-Material,  Heat  Conduction  Grid 
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Figure  6.  105-mm  XM797  Temperature  Recovery  Factor 


Figure  7.  105-mm  XM797  Heat  Transfer  Coefficient,  V=1.49  km/; 
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Figure  8.  105-mm  XM797  Heat  Transfer  Coefficient,  V=1.27  km/ 


Figure  9.  105-mm  XM797  Propellant  Temperature  Response 


Figure  10.  120-mm  XM797  Temperature  Recovery  Factor 
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Figure  11.  120-mm  XM797  Heat  Transfer  Coefficient,  V=1.65  km/s 
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Figure  12.  120-mm  XM797  Heat  Transfer  Coefficient,  V=1.43  km/s 


Figure  13.  120-mm  XM797  Propellant  Temperature  Response 


Contours  in  50  K  increments 


Figure  14.  120-mm  XM797  Temperature  Contours 
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Figure  15.  Effect  of  Chamber  Depth  on  120-mm  XM797  Propellant  Temperature  Response 
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APPENDIX: 

1-D  AXISYMMETRIC  HEAT  CONDUCTION  VALIDATION 
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This  appendix  presents  comparisons  between  numerical  and  analytical  results  for  a  1-D 
axisymmetric  heat  conduction  problem.  The  heat  conduction  code  described  in  the  main 
body  of  this  report  was  used  to  obtain  the  numerical  result.  The  configuration  is  an  infinite 
length  steel  cylinder  subjected  to  axisymmetric  convective  heating,  in  which  the  heat  transfer 
coefficient,  h,  is  held  constant.  The  problem  inputs  somewhat  typify  the  length  scale  and 
heating  environment  that  could  be  experienced  by  a  projectile  in  supersonic/hypersonic 
flight.  Table  A-1  shows  the  material  properties  of  the  cylinder,  as  well  as  the  initial  and 
boundary  conditions. 

Table  A-1.  1-D  Axisymmetric  Heat  Conduction  Validation  Problem  Conditions 


Material:  Mild  Steel 

Density,  p  (kg/m^) 

7850.0 

Conductivity,  k  (W/m/K) 

45.3 

Specific  Heat  Cp  (J/kg/K) 

444.6 

Diffusivity,  (m^/s) 

1.3x10-® 

Cylinder  Radius,  ro  (m) 

.0127 

Initial  Temperature,  Tj  (K) 

294. 

Adiabatic  Wall  Temperature,  Taw  (K) 

1800. 

Heat  Transfer  Coefficient,  h  (W/m^/K) 

5000. 

Nusselt  Number,  ^ 

1.4 

The  governing  equation  is 

d^T  1  ^ 

dr^  r  dr  a  dt 


An  analytical  solution  is  described  by  Schneider.^  The  temperature,  T,  is  a  function  of 
radial  location,  r,  and  time,  t,  and  is  given  by 


OO  1 


JljMn)  -Mie  -[  (M 

Mn  JiiMn)  +  mMn) 


(A-2) 


In  the  above  expression,  Jo  and  Ji  are  Bessel  functions  of  order  zero  and  one,  respectively. 
The  variable  0  is  defined  as 


(A-3) 


The  variable  M„  is  the  nth  root  of  the  eigenfunction  Bessel  equation 


M 


MM) 

MM) 


(A-4) 


*  Schneider,  P.,  Conduction  Heat  Transfer.  Addison- Wesley  Publishing  Co.,  Inc.,  Cambridge,  Massachusetts,  1955. 
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This  equation  has  an  infinite  number  of  roots,  which  are  found  using  Newton’s  method. 
The  equation  is  first  arranged  to  comprise  the  function  /: 


(A-5) 


The  derivative  of  /  is 


f'(M)  = 


M[Jl(M)  +  J?{M)] 
Jo 


(A-6) 


and  the  iteration  to  convergence  is  performed  using  the  algorithm 


=  Mi - 


mi) 


"  rm) 


(A-7) 


in  which  the  superscript  i  refers  to  the  iteration  number  for  any  particular  root  that  is  being 
sought.  The  initial  value  for  the  nth  root,  constructed  here  based  on  the  table  provided  by 
Schneider  for  the  eigenfunction  Bessel  equation  and  the  known  value  of  is 

Mi  =  1.4  for  n  =  1 

Ml  =  4.6  +  (n  -  2)ir  for  n  =  2, 3, 4, ...  (A-8) 


The  analytical  results  presented  here  were  generated  using  the  first  40  terms  from  Equation 
(A-2),  corresponding  to  the  first  40  roots  of  the  eigenfunction  Bessel  equation.  Forty  terms 
provided  accuracy  to  at  least  four  decimal  places. 

Figure  A-1  shows  the  temperature  response  of  the  cylinder  axis  (r/ro  =  0)  and  surface 
(r/ro  =  1)  as  a  function  of  time.  The  analytical  and  numerical  results  are  both  shown  and 
are  virtually  indistinguishable.  The  largest  observable  difference  between  the  numerical  and 
analytical  results  occurs  on  the  axis  at  approximately  1  second.  A  closer  examination  is 
made  in  Figure  A-2,  which  shows  the  numerical  and  analytical  temperature  distributions 
in  the  cylindrical  rod  at  t  =  0.9627  seconds.  The  agreement  on  the  axis  is  within  7  K, 
and  the  agreement  elsewhere  is  mostly  within  1  to  2  K.  The  disageement  near  the  cylinder 
axis,  which  is  acceptable  for  the  engineering  application  in  this  report,  is  attributable  to  the 
combination  of  large  grid  cells  near  the  axis  and  the  first  order  axis  boundary  condition  in 
the  numerical  code. 
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Figure  A-1.  1-D  Axisymmetric  Heat  Conduction  Validation:  Temperature  vs.  Time 


Figure  A-2.  1-D  Axisymmetric  Heat  Conduction  Validation:  Temperature  vs.  Location 
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LIST  OF  SYMBOLS 


Cp 

d 

h 

Jn 

k 

I 

Mn 

N. 

9 

r 

ro 

rj 

T 

t 

V 

X 


specific  heat 

propellant  chamber  depth 
surface  heat  transfer  coefficient 
Bessel  function  of  order  n 
thermal  conductivity 

reference  length,  equal  to  25.4  mm  (1  inch) 

nth  root  of  eigenfunction  Bessel  equation 

Nussclt  number 

surface  heat  transfer  rate 

radial  distance  from  x-axis 

radius  of  infinite  length  cylinder 

surface  temperature  recovery  factor 

temperature 

time 

projectile  velocity 

distance  along  projectile  axis,  measured  from  nosetip 


Greek  Symbols 

O'  thermal  diffusivity 

boundary  layer  transition  intermittency  factor 
p  density 


Subscript 

oo 

aw 

i 


I 

0 


t 


tr 

w 


free  stream  condition 
adiabatic  wall  condition 
initial  condition 

laminar  boundary  layer  flow  condition 
total  condition 

turbulent  boundary  layer  flow  condition 
transitional  boundary  layer  flow  condition 
wall  condition 
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US  ARMY  ARMAMENT  RD&E  CTR 
ATTN  AMSTA  AR  FSE 
E  ANDRICOPOULIS 
K  CHEUNG 
A  GRAF 
D  LADD 

PICATINNY  ARSENAL  NJ  07806-5000 

6  COMMANDER 

US  ARMY  ARMAMENT  RD&E  CTR 
ATTN  AMSTA  AR  CCL  B 
D  CONWAY 
D  DAVIS 
K  HAYES 
M  PINCAY 
F  PUZYCKI 
W  SCHUPP 

PICATINNY  ARSENAL  NJ  07806-5000 


3  U.S.  ARMY  RESEARCH  OFFICE 
ATTN  G  ANDERSON 

K  CLARK 
T  DOLIGOWSKI 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC  27709-2211 
2  DIRECTOR 

US  ARMY  BENET  LABORATORY 
ATTN  SMCAR  CCB  R 
P  AALTO 
S  SOPOK 

WATERVALIET  NY  12189 

1  DIRECTOR 

US  BELVOIR  RD&E  CTR 
ATTN  SATBE  FED 
N  BLACKWELL 
FORT  BELVOIR  VA  22060-5606 

2  COMMANDER 

UNITED  STATES  MILITARY  ACADEMY 
DEPARTMENT  OF  CIVIL 
AND  MECHANICAL  ENGINEERING 
ATTN  M  COSTELLO 
A  DULL 

WEST  POINT  NY  10996 

4  HERCULES  INC 

ALLEGANY  BALLISTICS  LABORATORY 
ATTN  J  CONDON 
W  NYGA 
J  PARRILL 
J  VILES 
PO  BOX  210 

ROCKET  CTR  WV  26726 


2  COMMANDER 

US  ARMY  ARMAMENT  RD&E  CTR 
ATTN  AMSTA  AR  CCH  A 
J  DIFUCCI 
M  PALATHINGAL 

PICATINNY  ARSENAL  NJ  07806-5000 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 


5  COMMANDER 

US  ARMY  ARMAMENT  RDLE  CTR 
ATTN  AMSTA  AR  FSF  T 
S  LIESKE 
J  MATTS 
F  MIRABELLE 
S  PUHALLA 
J  WHITESIDE 

41  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  SC  CC 
R  ANGELINI 
C  NIETUBICZ 
J  GROSH 
D  HISLEY 

AMSRL  SC  I W  STUREK 
AMSRL  WT  P 
A  HORST 
E  SCHMIDT 
AMSRL  WT  PA 
P  CONROY 
G  KELLER 
D  KRUCZYNSKI 
T  MINOR 
M  NUSCA 
AMSRL  WT  PB 
H  EDGE 
J  GARNER 
B  GUIDOS  (5) 

K  HEAVEY 
P  PLOSTINS 
J  SAHU 
P  WEINACHT 
AMSRL  WT  PC 
G  ADAMS 
R  FIFER 
AMSRL  WT  PD 
B  BURNS 
W  DRYSDALE 
R  LIEB 

S  WILKERSON 
AMSRL  WT  W 
C  MURPHY 
AMSRL  WT  WB 
F  BRANDON 
T  BROWN 
W  D’AMICO 
D  HEPNER 
M  HOLLIS 
AMSRL  WT  WC 
J  BORNSTEIN 
R  VON  WAHLDE 


AMSRL  WT  WD 
T  KOTTKE 
C  STUMPFEL 
AMSRL  WT  NC 
R  LOTTERO 
S  SHRAML 
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Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-MR-267 _ Date  of  Report  October  19Q5 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ _ _ .  _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ _ _ _ ' 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


